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The masses of higher D(nL) and D s (nL) excitations are shown to decrease due 
to the string contribution, originating from the rotation of the QCD string itself: it 
lowers the masses by 45 MeV for L = 2 (n = 1) and by 65 MeV for L = 3 (n = 1). 
An additional decrease ~ 100 MeV takes place if the current mass of the light 
(strange) quark is used in a relativistic model. For D S (1 3 D3) and D s {2Pi) the 
calculated masses agree with the experimental values for D s (2860) and D s (3040), 
and the masses of D(2 1 Sq), D(2 3 S\), D(1 3 D^), and D{1D2) are in agreement with 
the new BaBar data. For the yet undiscovered resonances we predict the masses 
M(D(2 3 P 2 )) = 2965 MeV, M{D(2 3 P )) = 2880 MeV, M(D(1 3 F 4 )) = 3030 MeV, 
and M(D S (1 3 F 2 )) = 3090 MeV. We show that for L = 2, 3 the states with j q = 1+1/2 
and j q = 1 — 1/2 (J = /) are almost completely unmixed ((f) ~ —1°), which implies 
that the mixing angles 6 between the states with 5 = 1 and 5 = (J = L) are 
9 40° for L = 2 and ss 42° for L = 3. 



I. INTRODUCTION 

Till recently only the low- lying IS and lPj states of the D and D s mesons were known 
from experiment [1]. The situation has changed recently owing to discoveries of new D s 
resonances: £> s (2710) [2, 3], £> s (2860) [2, 4], and D s (3040) [4]. Also last year, new D(L) 
states were observed by the BaBar Collaboration [5]: £>(2550), D*(2600), Dj(2750), and 
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-D}(2760), and in [6] the mass of D(1 3 P ) was measured with a good accuracy. 

The quantum numbers and decay modes of the new resonances were intensely discussed 
in a large number of recent studies [7]-[14] and also before in [15]-[25]. For £)* 1 (2710) the 
quantum numbers J p = 1~ were assigned [4, 7], although the analysis in Ref. [8] does not 
exclude that -D s i(2710) is an admixture of D s (2 3 Si) and D s (l 3 Di). The relatively narrow 
resonance D sJ (2860) with T = 48 ± 3 (stat) MeV is mostly considered as the L = 2 state 
with J p = 3~ [8]-[10], while in Ref. [25] the resonance with close values of the mass and 
width has the quantum numbers J p = + . The wide resonance _D s j (3040) is mostly assumed 
to be the 2P state with J p = 1+ [11]. 

The quantum numbers of the new D° resonances and their isotopic partners were dis- 
cussed in [12]-[14], where the broad resonance _D(2550) (T ~ 130 MeV) is considered as the 
singlet 2 1 S state, for which a large width as in experiment was obtained in Ref. [12], 
while much smaller total width was calculated in [13]. The resonance _D*(2600) with 
T = 93 ± 6 ± 13 MeV is consistent with the excited 2 3 Si state [14], or an admixture of 
the 2 3 5'i and 1 3 Di states with large mixing angle [12, 13]. The resonances _Dj(2760) and 
_Dj(2750) have relatively small widths, T ~ 60 — 70 MeV, and for them the quantum num- 
bers J p = 3~ and J p = 2~, respectively, were assigned in Ref. [13, 14], while in [19] (the 
second paper) these two resonances are considered as the same 1 3 Di state with J p = l~. 

These new data are extremely important for the theory to better understand the qQ 
dynamics and test predictions made in a large variety of models [15]-[25], some of which 
were made long ago [15, 16]. For low-lying states, the theoretical predictions are mostly 
in agreement with experiment within 20 — 50 MeV accuracy, although the parameters used 
may be very different. This is not surprising, because the very masses of low- lying states 
are usually used as a fit to determine the quark masses and parameters of the potentials. 
On the contrary, for higher states different predictions for the masses and the fine-structure 
(FS) splittings were obtained in different models. 

The new experimental data on the hyperfine (HF) splittings show that their values, 
~ 70 MeV, coincide for £> s *(2710) and £> s (2638) [26], £>*(2610) and £>(2540). The latter HF 
splitting was predicted in Ref. [27], if a "universal" coupling a H F = 0.31 is used in the HF 
potential. 

Also the experimental mass differences between £)*(2710) and D*(2610), D S (2Q38) and 
£>(2540), £> s (2860) and £>(2760), appear to be ~ 100 MeV, the only exception being the 
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D S (1P) multiplet, where on the one hand, the masses of D* 2 (2573) and 1^(2460), £> s i(2535) 
and A (2422) also differ by ~ f fO MeV, while M(D* s0 (2317) - M(D (2300) « 25 MeV and 
M(D s i(2460) - Di(2430) « 30 MeV are much smaller. Such small FS splittings cannot be 
explained within the "universal" description of fine structure used in Refs. [15, 16, 18]. 

This discrepancy has stimulated a lot of studies to understand why D* and D s \ have 
such small widths, Y < 3.8 MeV [1], and large mass shifts. Later it was understood that 
the dynamics of the D sJ (lP) multiplet is different for the states with the total angular 
momentum of the s quark j s = 1/2 and j s = 3/2 [21] -[25] and the bispinor structure of the 
D S (1P) wave function (w.f.) and the w.f. of the D(1S) mesons in the decay channel are very 
important. Two factors provide a large hadronic shift: the nearby S'-wave threshold and the 
large overlap integral between the upper components of the D s j(lP) w.f. with j s = 1/2 and 
the lower components of the w.f. of the D-meson in the decay channel [23]. 

Surprisingly, there are no large mass shifts for the other excited D and D s states observed 
and, as a whole, the single-channel description turns out to be a useful tool to understand 
the general structure of the D and D s spectra and FS splittings, and to predict the masses of 
the yet undiscovered resonances. Till now one of the best predictions for the meson masses 
of the low-lying states were obtained in the QCD motivated relativistic quark model (RQM), 
already in 1985 [15]. 

In contrast to the low-lying states, discrepancies show up for higher states, which may 
reach ~ 100 MeV. For example, for D s (3 1 S ) the masses 3097 MeV from the paper CTLS 
[7] and 3259 MeV [16], and for M(D S (2 S P 2 ) the values 3041 MeV [7], and 3157 MeV [16, 17] 
were obtained, showing differences > 100 MeV. The reasons why they occur will be discussed 
in the present paper. 

A comparison of the results obtained in different models is simplified, if the same value 
of the string tension a is used. The choice of a is of great importance, because the meson 
mass is proportional to y/a for the linear potential or, which dominates for higher states. 
However, much different values a are used in potential models: a large a ~ 0.26 GeV 2 in 
[16, 17] and small a = 0.115 GeV 2 in [24], a = 0.14 GeV 2 in [20]. Here we use a = 0.18 GeV 2 , 
which follows from the analysis of the Regge trajectories for light mesons, and was already 
used in Refs. [15, 19]. Taking the same a, one can establish common features and differences 
between the relativistic string Hamiltonian (RSH) [28] used here and the RQM developed 
in Refs. [15, 19]. In particular, we show that the choice of the current light (strange) quark 
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mass is of special importance in relativistic models. 

The only uncertainty in our calculations comes from the gluon-exchange (GE) potential, 
since at present there is no consensus about the value of the vector coupling at large distances, 
called the freezing constant or the critical constant a CT it- In Ref. [15] the value a cr it = 0.60 
was used and the variation of a cr i t in the range 0.60 ± 0.10 produces rather small changes, 
< 20 MeV, in the spin-averaged masses for higher states. However, the value of the strong 
coupling in spin-orbit and tensor potentials, ctps (/-*), is also not fixed now, in contrast to 
the FS in heavy quarkonia, where the scale \i and second order perturbative corrections are 
known [29], giving c*fs(/-0 smaller than a crit [30]. In our analysis of the FS here, we shall 
test different values of a^s to fit new experimental data on the D(1D) multiplet. 

In our paper we concentrate on the multiplets with L = 2, 3; to calculate mixing angles 
for states with J = L we use the basis j 2 from Ref. [22], where j g — I + s q is the total 
angular momentum of the light (strange) quark and the total spin J = j q + Sq is the sum of 
the light-quark total angular momentum and the spin of the heavy quark. It appears that 
for higher states with J = / the states with j q = 1 + 1/2 and j q = 1 — 1/2 are in fact unmixed, 
\4>{nl)\ ~ 1°, and this result may be important to study different decay modes of heavy-light 
mesons. Owing to the known relation between the mixing angle 4>{nl) and the mixing angle 
6{nL) in the S 2 scheme (or LS scheme with J = L + S), the states with the spin 5 = 1 
and S = (J = L) appear to be mixed with large mixing angle, e.g. 9 {ID) = 40.2°. 

II. RELATIVISTIC STRING CORRECTIONS 

Here we use the RSH, derived for spinless quarks and antiquarks [28], while all spin- 
dependent interactions are considered as a perturbation. To calculate the spectra of the 
heavy-light mesons this approach has some advantages as compared to the use of the Dirac 
equation (DE) and considering the heavy quark contribution as 1/mq corrections [16, 17]. As 
shown in Ref. [31], for a scalar potential the solutions of the DE have an important property: 
the spectrum is symmetric under the reflection of the eigenvalues (e.v.), e n — » — e n , so that 
negative energy states are in fact not present in the spectrum of a heavy-light meson. 

Moreover, from the expression for the squared e.v. of the DE (with a given / = l q and 
j = j q - the total angular momentum of a light quark) it follows that the mass difference 
between neighbouring states is equal to e 2 n+1 — e 2 n ~ Aa D + ln(^4M — ln(^K) [31], where 
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<7d is the string tension used in the DE and the constant k enters the Coulomb interaction, 
For the DE this mass difference (for a given a) appears to be significantly smaller than 
that in the RSH and the RQM, where it is equal to Ana. 

Just to compensate such a small spacing between radial excitations the larger value of 
the string tension, <td — 0.26 GeV 2 , is needed [16, 17] (in both cases the 1/rriQ corrections 
were taken into account). However, it remains unclear why in Ref. [16] the calculated values 
M{D S {\ 3 P )) = 2487 MeV and M{D s {lPf)) = 2605 MeV (for j q = 1/2) are similar to the 
numbers obtained in the RQM [15, 23], while much smaller values, 2325 MeV and 2467 MeV, 
were calculated within a similar approach in Ref. [17]. Here we will mostly compare our 
results with those models [15, 19], where the same o = 0.18 GeV 2 was used, and draw 
definite conclusions about the dynamics of the qQ interaction. 

The RSH H = H + H stT for spinless quarks and antiquarks was derived in instantaneous 
approximation [28] and has the following characteristic features: 

1. The QCD string, besides a standard rotation of a quark and an antiquark, rotates itself, 
giving an additional contribution to a Hamiltonian, H stT . For heavy-light mesons such 
string corrections are not large, ~ 30 — 70 MeV (for L — 1, 2, 3), and can be considered 
as a perturbation [18], while in light mesons the string corrections may dominate for 
states with large L [32]. 

2. In the unperturbed Hamiltonian H = T + Vq the kinetic term T [28] is 

T= ui ml uj 2 ml p 2 

2 2wi 2 2u 2 2w rcd ' 1 ' 

where by derivation the quark mass cannot be chosen arbitrarily and must be equal 
to the current mass rh q for the u, d, and s quarks and the pole mass mg for a heavy 
quark, thus taking into account perturbative corrections to the heavy quark mass. In 
our calculations m q = for the u, d quarks, m s ~ rh s (l GeV) = 200 MeV for the s 
quark, and the conventional pole mass m c (pole) = 1.42 GeV for the c quark [1] is used. 
This choice of m s is similar to that in Ref. [16], where m s = 220 MeV was used in the 
DE, being larger compared to m s {2 GeV) = 95 ± 20 MeV at the scale /i = 2 GeV [1]. 
The reason for that difference possibly originates from the fact that in the Hamiltonian 
approach the s-quark current mass m s {ji) enters at a smaller scale, fi ~ 1 GeV [33]. 
The value we take here, m s = 200 MeV, is significantly smaller than m s ~ 500 MeV 
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used in constituent quark models [19, 20]. It is important that the use of current quark 
masses allows to avoid several fitting parameters (constituent masses). 

3. The value of the string tension a — 0.18 GeV 2 cannot be used as a fitting parameter, 
as it is fixed by the slope of the Regge trajectories for light mesons. 

4. The choice of the GE potential is important for low-lying states. Here we use the 
vector strong coupling a B (r) which possesses the asymptotic freedom (AF) property 
and freezes at large distances at the value ct cr it- For higher excitations the choice of a cr it 
becomes less important; moreover, in many cases the GE potential can be considered 
as a perturbation. 

In the RSH H = H + H str the unperturbed part 

H = T(u 1 ,u 2 ) + V B (r) (2) 

contains the kinetic term T (1), where the variables have to be determined from the 
extremum conditions: ^ a = (i — 1,2) [28, 34]. Then one finds 

Ui(nL) = (y/p> + mt) nL (i = l,2). (3) 

The kinetic energy of a light (strange) quark is denoted as U\{nL) = u q (nL), and U2{nL) = 
u c is the kinetic energy of the c quark; the quantity cu rc d = J^!j^ 2 , and L — li + l 2 . Then 
putting Ui into Eq. (1), one arrives at a different form of T, denoted below as Tr: 

Tr = ^p 2 + m 2 q + y/p 2 + m 2 c . (4) 

Rigorously, the expression (4) for T R is valid only for L = 0, while in general, for L ^ 0, 
T = TR + T str contains the kinetic energy of the string rotation T str . For L < 4 this term T str 
is small compared to Tr and can be considered as a perturbation; its matrix element (m.e.) 
A str (nL) = (T stT ) nL is included in the mass formula (6). The form T R of the kinetic energy 
was suggested in Ref. [35] and used in many models [15, 36], while due to our derivation of 
T R one can establish the connection between the unperturbed RSH H and the RQM, where 
the same kinetic term is used. 

Then the e.v. M (nL) and w.f. are defined by the spinless Salpeter equation (SSE): 

[Tr + V B (r)\ VnL = M (nl)<p nL . (5) 
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TABLE I: The D and D s masses for the ID, 2P, and IF states (in MeV) 



StcltG 


GI[15] 


this 


\-J lAi \-J 


GXp. 


Refs [2-51 


£>(i 3 £> 3 ) 


2830 




2760 




2762 


D S (1 3 D 3 ) 


2920 




2840 




2860 


D S {2P?) 


absent 




3040 




3044 


D(1 3 F 4 ) 


3110 




3030 




absent 



It is essential that in the RSH approach the spin-averaged meson mass M(nL) = M cog (nL) is 
not only defined by the e.v. M (nl) (5), but also contains two additional negative contribu- 
tions: the string correction A str (nL) = {H stT ) nL [18, 34] and the nonperturbative self-energy 
(SE) term A SE (nL) [37]: 

M{nl) = M (nL) + A stI (nL) + A SE (nL). (6) 

For a given radial quantum number n, the string correction increases for larger L, while 
for a given L it decreases for higher radial excitations. For the D(1P), D(1D), and D(1F) 
states their values are equal to ~ —23 MeV, —45 MeV, and —65 MeV, respectively, which 
can be obtained using the analytical expressions for A str from Ref. [18]. As an illustration 
in Table I the masses calculated here for several D(nL) and D s (nL) states are compared 
to those from Ref. [15]. It appears that differences between them are mostly due to string 
corrections, ~ 40 MeV, and our numbers are closer to the experimental data [2]- [5]. 

For the IP and 2P states the string corrections are smaller, —22 MeV and —10 MeV, 
respectively, and the masses M(D S (2P 1 H )) calculated here coincide with the experimental 
mass of _D S (3040), if this resonance with J p = 1 + is identified as the higher 2P X H state 
with ] s = 1/2 (which has to have a larger total width), while in the low-mass state with 
M{2Pl) = 3020 MeV the state with j s = 3/2 dominates. 

Even larger mass differences occur for the yet unobserved states D(2 3 P 2 ), and 
D(n z F/±) (n = 1,2), for which the masses we predict here are ~ 100 MeV smaller than 
in [19] (see Tables II and III). 

The perturbative self-energy correction contributes to the current mass of a heavy quark 
(it gives ~ 15% for a c quark [1]) and moreover there exists a nonperturbative SE correction 
to the quark (antiquark) mass. This correction is very important to provide the linear 
behavior of the Regge trajectories [34]. As shown in Ref. [37], this correction is flavor- 



8 



dependent and strongly depends on the current quark mass, being small for a heavy quark 
and large for a light (strange) quark [37]: 

A SE = -—( Vf - VQ \ . (7) 
2n \Ld q (nL) ujQ{nL) J 

The factor i]f is determined by the quark current mass and the vacuum correlation length 
[37, 38]: r)f = 1.0 for a light quark, i] s = 0.70 for the s quark, and r\ c = 0.35 for the c quark. 
Notice that the number 3/2 enters the SE term (7), instead of the number 2 in [37]; this 
change follows from a more exact definition of the vacuum correlation length [38]. 

From Eq. (7) one can see that the kinetic energies Ui play a special role: they determine 
both the string and the SE contributions, and also enter all spin-dependent potentials [39] . 
In some potential models a negative overall constant Co is introduced, which may play the 
role of a self-energy correction, however, such a constant violates the linear behavior of the 
Regge trajectories; it is also important that in the RSH the SE terms decrease for higher 
states, being proportional to w~ 1 (nL). 

We use here the "linear+GE" static potential, Vb(t), which was already tested in a 
number of our previous works devoted to heavy-light mesons [18, 23] and heavy-quarkonia 
[40]: 

. . Act-sir) 
V B (r) =ar- (8) 

where the vector coupling o:b(V) is taken as in background perturbation theory [41] with 
ttcrit = 0.50, which is a bit smaller than a cr i t = 0.60 in Ref. [15], while a larger value 
«crit = 0.84 was used in Ref. [19]. In all cases the AF behavior of the vector coupling is 
taken into account. Notice that if a constant value a (without AF behavior) is used in the 
GE potential, then the value of a turns out ~ 30% smaller than a cr it- 



III. HIGHER D MESONS 



The masses of higher D excitations are presented in Tables II, III together with results 
from [15, 19], where the same a = 0.18 GeV 2 is used. In these Tables we have omitted 
results for the ground states, 1 1 S , l 3 Si, and the IP states, since they were studied in 
detail within the same approach in Ref. [18]; also for low- lying states their masses do not 
differ much in different models, since they are often used as a fit. 



9 



On the contrary, for higher states, large effect takes place when the constituent quark 
masses, instead of the current masses, are used. In Refs. [15, 19] the following masses were 
taken: 

Ref. [15] m uA = 220 MeV, m s = 419 MeV, m c = 1628 MeV, 
Ref. [19] m u4 = 330 MeV, m s = 500 MeV, m c = 1550 MeV, 
this paper, m u4 = 0, m s = 200 MeV, m c = 1420 MeV. (9) 

Our results are presented in Tables II and III. The mass M(D(1 3 D 3 )) = 2760 MeV cal- 
culated here, coincide with the experimental mass of £)(2760), which is assumed now to 
be the J p = 3~ state [13, 14]. This value is smaller than the masses 2863 MeV given in 
Ref. [19] and 2830 MeV given in Ref. [15] and this difference is partly explained by the string 
correction, equal to —45 MeV. 

Much larger differences occur for the excitations with n = 2 and L = 2,3. For example, 
M(2 3 L> 3 ) = 3212 MeV is obtained here, while the value 3335 MeV was predicted in Ref. [19], 
and this result cannot be explained by a string correction, which is only ~ —25 MeV in this 
case. From our point of view it happens due to the use of large constituent mass for a light 
quark. 

In Tables II and III we denote by P^ and Df the high-mass states with J = L, and 
by Pi and D\ the low-mass states. Each of these states is an admixture of the state with 
j q — I + 1/2 and j q — I — 1/2 in the j q 2 scheme, and for / = 2, 3 the mixing angle between 
these states appears to be small (see section V). 

For the 2 3 P2 state we predict the mass, 2965 MeV, smaller then the values 3012 MeV in 
Ref. [19] and 3035 MeV in Ref. [16]. For the states with L = 3 calculated here, the mass 
M(£>(1 3 F 4 ) = 3030 MeV, is 157 MeV and 80 MeV smaller than in Refs. [19] and Ref. [15], 
respectively, and these large differences can be only partly explained by the string correction, 
equal to —65 MeV for the IF states. The largest difference takes place here for M(2 3 F 4 ) = 
3430 MeV, which is much smaller than the value 3610 MeV from Ref. [19]. Again, such a 
large discrepancy cannot be explained by a string correction, which is ~ —48 MeV for the 
IF states. 

In Tables I-IV all FS splittings given are calculated taking the strong coupling a fs in the 
spin-orbit and tensor potentials equal to 0.45. 

We can summarize our results for the higher D mesons: 
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TABLE II: The D meson masses M(nS) (n = 2,3) and M(nD) (n = 1,2) (in MeV). The experi- 
mental data from [5]; the quark masses are given in (9); a = 0.18 GeV 2 in all cases. 



state 


exn [51 


this na.ner 


GI [15] 


EFG [191 


& lJ{) 


2539 


2567 


2580 


2581 


2 3 S 1 


2608 


2639 


2640 


2632 


3 1 .9n 


uiUOvll u 


3065 




3062 


3 3 5i 


CXi UOvll u 


3125 




3096 


l 3 Di 


a\ 
) 


2790 


2820 


2788 


1 3 D 3 


2763 


2760 


2830 


2863 


l£f 


absent 


2810 




2850 


1D 2 L 


2750 


2746 




2806 


2 3 L>i 


absent 


3215 




3228 


2 3 J D 3 




3212 




3335 



a ) The identification of this state is not certain, because the quantum numbers of the D(2760) 
state are not established. It could be either a 1 or a 1 3 Z?i state. 

TABLE III: The D meson masses M(nP),M(nF) (n = 1, 2) (in MeV) 



state 


this paper 


GI 


EFG 


2 3 P 

2if 

2P\ 

2 3 P 2 

1 3 F 4 

1 3 F 2 

2 3 F 4 


2880 
2960 
2940 
2965 
3030 
3088 
3430 


3110 


2919 
3021 
2932 
3012 
3187 
3090 
3610 



1. The HF splitting between D(2 3 Si) and D(2 1 So), equal to 72 MeV, was calculated 
with the use of the "universal " strong coupling in the HF potential, q;hf — 0.31 from 
Ref. [27]; this splitting is in full agreement with experiment. 

2. The string corrections, present in the RSH, reduce the spin- averaged masses of the nL 
multiplets by -25 MeV, -45 MeV, -65 MeV for the states with L = 1, 2, 3, respectively. 
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3. Large mass differences for high excitations like D(2 3 D 3 ) and D(2 3 F 4 ) reach 120 MeV 
and 180 MeV compare to the predictions in Ref. [19]. 

4. The recently observed D(27Q0) and _D(2750) resonances are interpreted as the D(l 3 _D 3 ) 
and the low-mass D(1D\) states, where DQDY) is in fact the state with j q = I + 1/2 
(see section V) and therefore should have relatively small total width, as it is observed 
in the BaBar experiments [5]. It implies that in the LS scheme the states D{1 3 D 2 ) 
and D(l 1 D 2 ) are mixed with the mixing angle 8 = 40°. 

IV. HIGHER D s MESONS 

For the S-wave states, there are no string corrections, nevertheless, the mass 
M(D S (3 1 S ) = 3140 MeV calculated here, is 79 MeV less than the one given in Ref. [19] (see 
Table IV). From our point of view, this happens because of the large constituent mass m s = 
500 MeV taken in Ref. [19]. To illustrate this effect we have solved the SSE with two different 
masses of the s quark: m s = 0.5 GeV and m s = 0.2 GeV, keeping all other parameters the 
same. Then the mass difference 8(nL) = M cog (nL, mi = 0.5 GeV)— M cog (nL, mi = 0.2 GeV) 
appears to be almost constant for a fixed L and changing n: 5(2P) ~ 5(3P) = —138 MeV; 
5(1D) ~ S(2D) = -130 MeV, and 5 (IF) ~ S(2F) = -120 MeV, and thus one may expect 
mass differences ~ 100 — 150 MeV to occur between relativistic models with large constituent 
light (strange) quark mass compared to the RSH, which uses curent-quark masses. 

Just for that reason the masses M(D S (2 3 D 3 ) and M(D S (1 3 F 4 )) are in our calculations 
~ 180 MeV lower than in Ref. [19] and again such a large difference cannot be explained by 
the string corrections, which is only ~ —45 MeV for the D S (2F) state. 

Thus one can conclude that a large decrease in the masses of higher states predicted here, 
mainly comes from two sources: the string correction and the use of the current mass for an 
s quark, which is significantly smaller than a typical constituent mass m s ~ 450 ± 50 MeV. 

There exists another characteristic feature of the D s spectrum - for all known states the 
experimental masses of D s (nL) and D(nL) differ by ~ 100 MeV. In our calculations such a 
spacing 8 s (nL) comes from two sources; first, from different e.v. M (nL) of Eq. (5) in the 
cases with m q = and m s = 0.20 GeV, which gives ~ 50 ± 10 MeV difference. Secondly, 
the light and the s quarks have different nonperturbative SE corrections (negative), which is 
~ 40 ± 10 MeV smaller for the s quark as compared to a light quark. Altogether S s appears 
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TABLE IV: The masses M(nL) (in MeV) for D s mesons 





exD [11- [41 


this 


l_/ l_/ C 


GI [15] 


EFG [191 


2 1 Sn 


2638 a 




2656 


2670 


2688 


2 3 1 S'i 

— ■ kJ \ 


271(1 

2688 b 




2728 


2730 


2731 

1 til 


3 1 .S'n 

O A_>(J 


Cki UOC11U 




3140 




3219 




oKopvif 
C I UOCllt 




3200 




3242 


2 3 Pn 


oKopyif 




2970 




3054 


2P, H 


3044 




3040 




3154 


2P, L 






3020 




3067 


2 3 Pi 






3045 




3142 


1 3 £)i 


oKopvif 




2870 


2900 


2913 


ID? 






2885 




2961 




oKopvif 




2828 




2931 


1 3 D? 


2862 




2840 


2920 


2973 


2 3 Di 






3290 




3383 


2 3 A3 


absent 




3285 




3469 


1 3 F 4 


absent 




3110 


3190 


3300 


1 3 F 2 


absent 




3150 




3230 


2 3 F 4 


absent 




3490 




3754 



The data of SELEX [26]. b The data of Belle [3]. 



to be ~ 100 MeV for low-lying states and smaller, 5 S ~ 70 — 80 MeV, for higher states. 
Our results about the D s spectrum can be summarized as follows 

1. The HF splitting between D s (2 3 Si) and D S (2 1 S ), calculated with the use of the 
"universal" a^F = 0.31 from Ref. [27], gives good agreement with experiment. 

2. The resonance _D S (3044) is considered here as the high-mass state D S (2P 1 H ), which is 
dominantly the state with j s = 1/2 (see section V) and therefore has to have large 
total width, in agreement with experimental value T = 239 ± 35 MeV; also its mass, 
3040 MeV, calculated here, is in full agreement with experiment. 
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3. The resonance _D S (2860) is interpreted as the D S (1 3 D 3 ) state and its calculated mass 
2840 MeV is in agreement with experiment. This state with J p = 3~ and j s = 5/2 is 
assumed to have relatively small total width, as it takes place for D* 2 (2573). Indeed, 
the experimental width r(D(2860)) = 48±3 MeV [2, 4] is small for so high a resonance. 

4. The calculated masses of the higher states, like D S {2D) and D S {2F) are 120-200 MeV 
less than those from Ref. [19]. 

V. FINE STRUCTURE SPLITTINGS 

On a fundamental level, the spin-dependent (SD) potentials K(r) (i — 1 — 4) have been 
studied in analytical approaches [39, 42], and also on the lattice [43], where the SD potentials 
are expressed via the vacuum correlators. When the spin-orbit potential Vso(r) is considered, 
its perturbative part can be expressed only through the vector potential V 2 (r) = V(r) if 
the Gromes relation [44] is used, and its nonperturbative part is determined by the scalar 
confining potential S(r) = or. For the tensor potential the nonperturbative contribution 
appears to be very small [39, 43] and it is defined by the perturbative potential only, usually 
denoted by V 3 (r), which in general case is not equal to [^ — V"~\ , as it takes place for the 
one-gluon-exchange (OGE) potential (notice, that in the static potential (8) the effective 
vector coupling a^ir) includes higher order perturbative corrections, while these corrections 
appear to be different for different spin-dependent potentials and in OGE approximation 
they are neglected): 

V s (r) = 3T (r)£ = t (10) 

Here the factor £(nL) is introduced to show the difference between V 3 (r) and 3T (r). In 
heavy quarkonia (HQ) this factor £ appears due to second order perturbative corrections, 
being £ ~ 1.30 ± 0.05, both in charmonium and bottomonium [30]. However, the value 
of £ remains unknown for heavy-light mesons and the difference between Vs(r) and 3To(r) 
may be important for the FS analysis. In HQ for the IP states the spin-orbit ago(l-P) and 
the tensor t(lP) m.e. can be extracted from the experimental masses since all of them are 
measured with great accuracy. 

The study of the FS of the D{nL) and D s {nL) multiplets is a more complicated task, 
because only a few masses are known from experiment, besides those for the D(1P) and 
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D S (1P) multiplets. Moreover, many states lie above open thresholds and may have mass 
shifts, which change the mass values as compared to those in single-channel approximation. 
Nevertheless, a general analysis of the FS in heavy-light mesons is very useful and allows to 
understand better the FS dynamics and make definite conclusions about mixing angles for 
the states with J = L. 

For a multiplet nL the FS is considered here in the basis j 2 , where the total angular 
momentum of a light (strange) quark j q is diagonal [22]. (Below we use the notation j q = j). 
This basis is especially convenient for the calculation of the mixing angle (denoted as 4>{nl)) 
between the states with j = 1 + 1/2 and j = I — 1/2, if J = I. Another scheme, S 2 , is also 
often used, and in this scheme the notation 9 for a mixing angle is used here. The relation 
between 9 and can be easily established, writing the high-mass state (Vj) and low-mass 
state (Lj) with J = L in both schemes. In the j 2 basis we write 

|J H ) = sin0|j = / + ^) + cos0|j = /-i), 

|J L ) = C os0|j = / + i)-sin0|j=/-i), (11) 
while in the S 2 scheme the same physical states are defined as in Ref. [12], 

|L H ) = -sm9\ 1 Lj) + cos9\ 3 Lj), 

\L h ) = cos 01%) +sinfl| 3 Lj). (12) 
Then taking from Ref. [22] the relations: 



|J = /,j = /-l/2) = \/^l^ = ^=i>-V^I^ = ^ = 0>> 



|j = /,j = / + i/2) = V^TTl J = / ' 5 = 1 ^ + V27TT |J = / ' 5 = 0) ' (13) 

and inserting them into Eq. (11), one obtains 



9 = -0 + arccos-y/^^. (14) 

For L = 1,2,3 it gives 

9(L = 1) = -0 + 35.26°, 9(L = 2) = -0 + 39.23°, 9(L = 3) = -0 + 40.89°. (15) 

To determine one needs to know the m.e. of the spin-orbit and tensor potentials, which 
are written here in a more general form than in Ref. [22]: 

Vso = A 1 (r)2rs 1 + A 2 (r)2Z-s 2 , (16) 
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with 

w , 1 V'-S' IV w , 1 V'-S' IV 

Al(r) = 4^— +2^- A2(r) = 4^f^— + 2^T- (17) 

Notice that the kinetic energies Ui enter Aj in Eq. (17) instead of the constituent masses 
usually used in potential models. This change follows from the general consideration of 
spin-dependent potentials in the RSH [39] and is important for higher states, decreasing 
their FS splittings. 

For the linear confining potential S' = a, while the perturbative vector potential V(r) is 
taken in the form, satisfying the relation V'(r)/r = 4aps/r 3 = T , as for the OGE potential, 
where the vector coupling aps is considered as an effective coupling. Then the quantity V3 
in the tensor potential, 

is given in Eq. (10) and the tensor operator is defined as usual by 

= 3 ( g ,TK^T) _ ^ (lg) 



Later we use for simplicity the notations \i(nl) for m.e. (Aj(r)) n j and 

\3wiw 2 / ni 



(20) 



The spin-orbit m.e. aso, given by 

a so (nl) = Ai(raZ) + A 2 (nZ), (21) 

and the tensor m.e. t(nZ) fully determine the FS splittings for the states with J — I + 1 and 
J = I — 1 (in both cases spin 5 = 1): 

M(J = l + l,S = l) = M cog + la so - 2(2 / +3) *» 

M( J = Z - 1, 5 = 1) = M cog - (Z + l)a so - 2(2f + _\ ) t ( 22 ) 

For J = Z the states with j = 1 + 1/2 and j = 1 — 1/2 are mixed and their masses and mixing 
angles are defined by the matrix M mix : 

M . = I + ^fey [t - 8A 2 (Z + 1)] -(4A 2 - j ^ 

-(4*2 -f)^ -a5o(/ + l) + 5^iy(« + 8ZA 2 ). ' 
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From Eq. (23) one can see that in general the matrix M mix depends on a so and t, and also 
on the m.e. A2, and the value of the factor 4\ 2 — t, present in the off-diagonal m.e., is 
important for the determination of the mixing effect. 

In the heavy-quark limit there is no mixing, because both A2 and t are going to zero 
(they are proportional to ra^ 1 (n = 1,2) and may therefore be neglected). Then the high- 
mass state H has j — I + 1/2 while the low- mass state has j — I — 1/2, if a so is positive. 
However, such a situation with a S o ^ t does not occur even in bottomonium, where for the 
IP states a so (bb,exp) = 13.65 ± 0.39 MeV coincides with t(66,exp) = 13.13 ± 1.04 MeV 
within the experimental errors, and in charmonium aso(cc, exp) = 34.96 ±0.13 MeV is even 
14% smaller than t(cc,exp) = 40.63 ± 0.26 MeV (their ratio is 0.86). 

Such a decrease of the spin-orbit m.e. and the ratio a so /t occurs due to the negative 
—<j{r~ l ) nL term and a partial or full cancellation in the m.e. ( v '~ s ' ) n L is possible. (Also 
the m.e. (r~ 3 )„£, entering the spin-orbit and tensor m.e., decreases for increasing n and L). 
Therefore it is of interest to define the quantity 

4 

A so (nL) = -a FS (r- 3 ) nL - <7<r _1 ) n L, (24) 

which does not depend on and enters \ = + ^ (i — 1,2). In the D(D S ) mesons 
the factor Aso(nL) is negative and its magnitude depends on the value of «fs taken. Here 
a FS = 0.45 is mostly used, which is a bit larger than a so ~ 0.38 ± 0.02 extracted from the 
charmonium FS [30]). The values of A so can be illustrated by the following numbers: 

1. In bottomonium Aso — 0.14 GeV 3 is positive and relatively large, while in charmonium 
A so = ±0.01 GeV 3 is already small, even compatible with zero, so that the ratio 

= 0.86 is less than unity. 

2. For the D(nL) multiplets the factor Aso is always negative: ~ —0.017 GeV 3 for the 
1P,2P states, and ~ -0.028 GeV 3 for the ID and IF states (for a FS = 0.45). This 
result weakly depends on the quark masses used. 

3. In M mix a common scale is defined by the tensor m.e. t, and for t = T (i.e., £ = 1.0) 
it has values: t(lP) = 39 MeV, t(2P) = 29 MeV, t(lD) = 11.3 MeV, t(lF) = 5 MeV. 
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4. For L = 2,3 the mixing angle is very small, \(f>\ < 1°, for any reasonable choice of 
coupling. On the contrary, for the nP (n = 1,2) states the mixing angle is very 
sensitive to a^s used. 

5. For small coupling, a F s < 0.30 the mixing angle decreases, so that the main un- 
certainty in any FS analysis comes from the value of aps taken, which is not fixed 
yet. 

In our calculations the following values of the kinetic energies are obtained for D(nL): 

u q {lP) = 0.60 GeV, u q (lD) = 0.683 GeV, u q (lF) = 0.757 GeV, 

u c (lP) = 1.555 GeV, u c (lD) = 1.588 GeV, u c (lF) = 1.62 GeV. (25) 

For the D s mesons the FS picture is essentially the same, because the m.e. for D s , which 
are important for the FS, coincide within 1-5% with those of the D mesons, and therefore 
the D s FS splittings and mixing angles are practically the same as for the D mesons (see 
Tables II-IV). 

We also assume that for a given nL multiplet the masses of the M(J — 1 + l,j — 1 + 1/2) 
and M(J = l,j = 1 + 1/2) states have no mass shifts (or have small mass shifts), as it 
happens for the D(1P) and D S (1P) multiplets, and therefore the mass differences between 
these states, 

M(D* s2 (2573) )- M(D sl (2535) = 37.31 ± 1.0 MeV, 
M(D* 2 (2A60) - M(£> 1 (2422)) = 40.8 ± 1.6 MeV, (26) 

may be considered as the most stable characteristic of a given multiplet nL; in general this 
mass difference is denoted by A(raZ): 

A(j = I + 1/2) = M(J = I + 1, j = I + 1/2) - M(J = I, j = I + 1/2). (27) 

Our calculations show that for the D(1P) states the quantity Ai = 5.5 MeV is positive 
and small, while A2 = 17.3 MeV is relatively large, giving aso(l-P) — 22.8 MeV, and 
t(lP) = T (1P) = 38.6 MeV (f = 1.0), so that the ratio a so /t = 0.59 is smaller than in 
charmonium, where this ratio is 0.86. 

With the use of M mix , Eq. (23), the mass splittings within the 2P multiplet are calculated 
(see Tables III and IV) and the mixing angle depends on a FS , decreasing for smaller coupling: 
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a large angle (f>(lP) = —38° is obtained for large a FS = 0.60, while <f>(lP) = —12.6° for 
a FS = 0.45 and = —4.2° for a smaller a FS = 0.33. 

For higher orbital excitations (1 — 2,3, n = 1) the nondiagonal terms in the matrix M m i x 
appear to be much smaller than the diagonal m.e. for ctpg = 0.45 and due to this fact the 
mixing angles 

<t>{lD) = -0.89°, 4>{1F) = -1.0° (28) 

are very small. From Eq. (14) these values of correspond to the following angles 9 between 
the states n 3 Lj (S = 1) and n x Lj (S = 0) with J = L: 6 (ID) = 40° and 6 (IF) = 42°. 
This result may be important for the hadronic decays of these resonances [12]. 
The calculated mass differences A(nL), defined in Eq. (27), 

A(LP) = 37.6 MeV, A(l£>) = 15 MeV, A(1F) = 14 MeV, (29) 

are in good agreement with the experimental numbers: A(1P, exp) = 37.6 MeV [1] and 
A(ID, exp) ~ M(L>(2760)) - M (D(2750)) = (11 ± 9) MeV [4]. 

We do not discuss here the masses of the states with j — I — 1/2, which may have large 
mass shifts. In the single-channel approximation the mass M(1P, + ) is 104 MeV smaller 
than M(1P, 2 + ); for L = 2 almost equal masses M(l 3 D 3 ) and M (1 3 Di) are obtained, while 
for the L = 3 states M(l 3 F 4 ) is even smaller than M( 3 F 2 ) (see Tables II, III). It does not 
exclude that because of possible mass shifts, the physical masses M(l 3 Di) and M(1 3 F2) 
become smaller than M(1 3 D 3 ) and M(1 3 F 4 ). 

In Ref. [22] for the IP states the approximation A 2 = § was used, which in our consid- 
eration is also valid for the IP and 2P states. For the ID and IF states the factor A 2 is 
smaller, X 2 (1D) ~ 0.3t, A 2 (l-F) ~ 0.2t, and therefore the factor 4A 2 — t in the off-diagonal 
term in Eq. (23) is also smaller. However, the main reason why a small mixing occurs for 
I = 2, 3, is that the diagonal terms appear to be larger than the off-diagonal terms due to 
larger algebraic coefficients. 

As a result, for / = 2 and 3 the high-mass state is dominantly the state with j — I — 1/2 
and the low- mass state is mostly the state with j — I + 1/2. 
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VI. CONCLUSIONS 

The spectra of the D and D s mesons were studied with the use of the RSH, where only 
such fundamental parameters as the string tension and the quark current masses enter, and 
the only uncertainty comes from the freezing constant a CT it, which for higher states gives a 
theoretical error < 20 MeV in the spin-averaged mass. We have shown that 

1. The calculated masses of the higher excitations appear to be 50-150 MeV lower than 
in other RQM with the same string tension a — 0.18 GeV 2 . It occurs for two reasons: 
first, due to the string corrections for the states with L ^ and secondly, because we 
use the current quark masses. 

2. Using the j 2 basis, the states with J = I are shown to have very small mixing angles 
for I — 2 and 3: <p ~ —1°. It means that the states 1 3 L^ and 1 1 Ll are mixed with 
9(1D) = 40° and 9 (IF) = 42°. 

3. The calculated masses of the state 1 3 D 3 and the low-mass state ID 2 agree with the 
new BaBar resonances, _D(2760) and £>(2750). 

4. The resonance D*j(2860) is considered as the 1D S (1 3 D 3 ) state and £> s (3040) as the 
high-mass 2P 1 H state. 

5. For the yet unobserved resonances the following masses are predicted: M(D(2 3 P2)) = 
2965 MeV, M(2P X L ) = 2940 MeV M(D(l 3 F A )) = 3030 MeV and M(D s (l 3 F A )) = 
3110 MeV. 
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